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ABSTRACT 



> 

_ Aims. Analyses of recent cosmic microwave background (CMB) observations have provided increasing hints that there are deviations in the 
I . universe from statistical isotropy on large scales. Given the far reaching consequences of such an anisotropy for our understanding of the 
' universe, it is important to employ alternative indicators in order to determine whether the reported anisotropy is cosmological in origin and, if 
so, extract information that may be helpful for identifying its causes. 

Methods. Here we propose a new directional indicator, based on separation histograms of pairs of pixels, which provides a measure of 
departure from statistical isotropy. The main advantage of this indicator is that it generates a sky map of large-scale anisotropies in the CMB 
^ ^ temperature map, thus allowing a possible additional window into their causes. 

I Results. Using this indicator, we find statistically significant excess of large-scale anisotropy at well over the 95% confidence level. This 

O anisotropy defines a preferred direction in the CMB data. We discuss the statistical significance of this direction compared to Monte Carlo data 

obtained under the statistical isotropy hypothesis, and also compare it with other such axes recently reported in the literature. In addition we 

^ show that our findings are robust with respect to the details of the method used, so long as the statistical noise is kept under control; and they 

. . remain unchanged compared to the foreground cleaning algorithms used in CMB maps. We also find that the application of our method to the 

^ first and three-year WMAP data produces similar results. 
• i-H 

d Key words. Cosmology - large-angle anisotropies - cosmic microwave background - large-scale structure of the universe 



1. Introduction 

Recent years have witnessed a tremendous accumulation of 
high-precision cosmological data. In particular, the combina- 
tion of the first-year and, more recently, the three-year high- 
precision data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) (Bennett et al. |2003a[ |2003bl Hinshaw et 
al. 120061 ) has produced a wealth of information about the 
early universe. These data have provided evidence of a nearly 
spatially-flat universe with a primordial spectrum of almost 
scale-invariant density perturbations, as predicted by generic 
inflationary models. 

On large scales, however, a number of potentially impor- 
tant anomalous features in the cosmic microwave background 
(CMB) have been reported (see Copi et al. |2005l |2006! for a 
detailed discussion), including their low quadrupole and oc- 
topole moments (Spergel et al. 120031 120061) . the quadrupole- 
octopole axis alignment (Tegmark et al. 120031 de Oliveira- 



Costa et al. 120031 Weeks 120041 Bielewicz et al. |2005l Wiaux 
et al. 120051 Abramo et al. 120061 ). evidence of a North-South 
asymmetry (Hansen et al. I2004al Eriksen et al. I2004al Bernui 
et al. ,2006) . indications of a preferred axis of symmetry around 
(b - 30°, I - 250°) in galactic coordinatesjj or directions of 
maximum asymmetry towards {b - 100°,/ = 237°) (Bunn 
& Sc ott 2QQQ, Copi et al.|2 004l E riksen et al. I2004at Hansen 
et al. I2004al I2004bl Wibig l2005l Land & Ma gueijo 
I2005bl Schwarz et al. [20041 Prunet, et al [20051 Donoghue & 
Donoghue 120051 ). as well as indications of non-Gaussian fea- 
tures in the CMB temperature fluctuations (Copi et al. |2004t 
Komatsu et al. ' 20031 Park [20041 Chiang et al. 20W. Vielva et 
al. ,2004. Eriksen et al. [20051 Larson & Wandelt,2004, Coles et 
al. [20041 Naselsky et al. HOQl Bernui et al. [2006b . 

Clearly the study of such anomalies must take into account 
that they may have non-cosmological origins such as unsub- 
tracted foreground contamination and/or systematics (Schwarz 



Send offprint requests to: M.J. Rebougas 



' Throughout this paper we use galactic coordinates with equator 
defined by = 90°, / € [0°, 360°]. 
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et al. 120041 ). unconsidered local effects like gravitational lens- 
ing (Vale 120051 see, however, Cooray & Seto |20051 l. or other 
mechanisms to break statistical isotropy (Gordon et al. 120051 
Tomita 12005). They may also have an extra-galactic origin 
(Bunn & Scott 2000; Tegmark et al. l2003l Eriksen et al. l2004al 
Hansen et al. 2004ai Land & Magueijo|2005i 12005 b'). If they 
turn out to have a cosmological origin, however, this could 
have far-reaching consequences on our understanding of the 
universe and, in particular, on the standard inflationary picture 
that predicts statistically isotropic CMB temperature fluctua- 
tion patterns and Gaussianity. 

In view of this, a great deal of effort has recently gone into 
verifying the existence of such anomalies by employing sev- 
eral different statistical approaches. Apart from coiToborating 
the existence of a large-scale anisotropy, using multiple indi- 
cators may be useful in determining their origins. In addition, 
different indicators can in principle provide information about 
the multiple types of anisotropy that may be present (Souradeep 
& Hajian 2004, 2005; Hajian & Souradeep|2003l|2005l Hajian 
et al. 120041 ). 

Here we propose a new indicator, based on the angular pair- 
separation histogram (PASH) method (Bernui & Villela 2006 ), 
as a measure of large-scale anisotropy. An important feature 
of this indicator is that it generates a sky map of large-scale 
anisotropics for a given CMB temperature fluctuation map. 
This level of directional detail may also provide a possible ad- 
ditional window into their causes. 

The structure of the paper is as follows. In Sect. |2] we in- 
troduce our anisotropy indicator. Section[3]contains the results 
of the applying our indicator to the WMAP data, and finally 
Sect. |4] contains the summary of our main results and conclu- 
sions. 



2. Anisotropy Indicator 

In this section we construct an indicator that could measure the 
departure of CMB temperature fluctuation patterns from sta- 
tistical isotropy. The construction of this indicator involves a 
number of steps. To begin, we recall that in CMB studies the 
celestial sphere is discretized into a set of equally sized pix- 
els, each with a temperature fluctuation value. The first step 
involves the ordering of the pixels by temperature. The pixels 
are then separated into two sets, one with positive and the other 
with negative temperature fluctuations, and each set is subdi- 
vided into a number of submaps, each consisting of an equal 
number of consecutive pixelso This ordering and subdivision 
of the CMB map is motivated by the fact that any anisotropy 
in the CMB temperature fluctuations implies the presence of 
uneven distribution of the pixels with similar temperatures. 

To quantify this distribution, the next step is to generate 
the pair angular-separation histograms (PASH's), which are ob- 
tained by counting the number of pairs of pixels in a given 
submap with angular separation a lying within small sub- 
intervals (bins), Ji, of (0, 180°], of length 6a = 1807A^bins' 
where 7,- - (a,- - ^ , -H , / = 1,2,..., N^^^^ , with the 

^ We shall see below that, subject to certain reasonable constraints, 
the details of this subdivision will not change our results significantly. 



bin centers at a, - (/ - j)6a . The PASH is then defined by the 
normalized function 



<!>(«,) 
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where n is the number of pixels in the submap, T](a) the num- 
ber of pairs of pixels with separation a, and where the normal- 

ization condition 2,=i ^ ^(ad 6a - 1 holds. In this way the 
pixels in each submap are ordered in a histogram. Having cal- 
culated the PASHs for the submaps, we average them for each 
bin to obtain the mean PASH (MPASH) ( <D(a,) ) . 

To obtain a measure of anisotropy for an observational pix- 
elized map, we compare the MPASH calculated in this way 
with the histogram expected from a statistically isotropic map. 
To obtain this map, we first use the expected number of pairs 
^exp(Q')^ with angular separation a e 7,-, for a distribution of 
n pixels in the sky sphere 5'^. This allows the normalized ex- 
pected pair angular-separation histogram (EPASH) to be writ- 
ten as 



<l>exp(Q',) = ^ ^ Z ''exp(«) = ^ ^exp(a,) , 



(2) 



where - n(n - l)/2 is the total number of pairs of pixels, 
'Pexpiad = YjaeJi Vexpio:)/N is the expected probability that a 
pair of objects can be separated by an angular distance that lies 
in the interval 7,, and where the coefficient of the summation 
is a normalization factor. Defined in this way, (i>expiai) gives 
the probability distribution for finding pairs of points (pixels 
in a submap) on the sky sphere with angular separations a, € 
(0, 180°]|1 

We denote the difference between the MPASH, 
( <l)Q|^g(a,) ), calculated from the pixelized observational 
map, and the EPASH <I>exp(ffi) , obtained from an statistically 
isotropic map, as 



T(ff,) = <3>obs(ff,))- Oexp(a/) . 



(3) 



The latter can be obtained from an analytical expression de- 
rived by Teixeira (2003 ), which gives the expected probability 
that an arbitrary pair of pixels, in a spherical cap of aperture yo 
degrees, is separated by an angle y [ y e (0, 2yo ^ 180° )] 



!^exp(r; ro) 



siny 



( arcsin[ {(cos yo + cos y) 



;r(l - cosyo)^ 
-^cosy - cos 2yo ) / {sin yo ( 1 + cos yo) + cosy ) ] 

+ (1 - 2 cosyo) arccos[cotyo tan(y/2)] ) . (4) 

Alternatively, one can approximate <l>exp by the MPASH, <<!>), 
obtained through the average of the PASH's produced from a 



^ In this connection it is instructive to note that, for a homogeneous 
distribution of points on 5^, all angular separations < y < 180° 
are allowed, and the corresponding probability distribution can be cal- 
culated to give ?'gxp ^^(y) = I siny. This represents the limit of 
a statistically isotropic distribution of points in 5^ as the number of 
points go to infinity. One can thus quantify departure from statisti- 
cal isotropy by calculating the departure of the mean observed prob- 
ability distribution ( OoftjC"^;) ) from this quantity, i.e. by evaluating 

< ^obMi) > - <t>expiai) . 
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Fig. 1. Plots of the difference histogram T{a) illustrating the 
magnitudes of the deviation from isotropy as a function of sepa- 
ration angle calculated for two specific spherical caps, centered 
at {b = 65°, I = 56°) and {b = 1 15°, / = 236°), with a yo = 30° 
of aperture. Plotted also, as a continuous line, is the mean T 
corresponding to the set of spherical caps into which the sphere 
was subdivided, along with the corresponding two-standard de- 
viation region as a shaded area. As can be seen, the measured 
anisotropy is significant in one cap {{b - 115°,Z = 236°)] but 
not in the other. 



set of Monte Carlo-generated statistically isotropic maps (see 
Bernui & Villela|2006]for details). 

Clearly, the difference histogram T gives the departure of 
the mean angular-separation probability distribution, obtained 
from observational data ( ^q^^ ), with respect to the probabil- 
ity distribution of the statistically isotropic case. In this way, 
it contains information about the amplitude and angular scale 
of the large-angle correlations for the sky sphere, but no di- 
rectional information about possible anisotropics. However, T 
can also be calculated for sections of the sphere, in which 
case it can be used to study the differences in the extent of 
anisotropy in different regions of the sky. As an illustration of 
this point, we have plotted in Fig. [T] difference histograms T 
corresponding to two different (antipodal) spherical caps cen- 
tered at {b = 65°, / = 56°) and {b = 1 15°, / = 236°), each with 
an aperture of yo = 30°. These histograms illustrate the magni- 
tude of the deviation from isotropy as a function of separation 
angle. As can be seen, the large-scale angular correlations are 
significant in one cap but not in the other, in comparison with 
the 2-standard deviations error bars. 

The study of the individual caps can thus provide hints of 
possible large-scale anisotropics. But for a more systematic 
study, one needs to understand how the anisotropy varies di- 
rectionally over the whole sphere. One therefore requires an 
indicator that encodes a measure of anisotropy as a scalar func- 
tion defined over the celestial sphere, rather than a histogram 
for each cap. 

To define this indicator, let = Q.(6j,<pj;jo) e he 
a spherical cap, with an aperture of yo degrees, centered at 
(0j, (pj), for j - 1, . . . , A^caps, where the union of the A^caps 
caps covers the whole celestial sphere S^. Define the scalar 
function o" : "R^, that assigns to the 7*-cap, centered 



at (6j, (fij), a real positive number crj given by the variance of 
the T (which has zero mean) 

%ns j-t 

where = T(Oj, (pj) is the T for the cap centered at {6j, (pj). 

To obtain a local directional indicator, we cover the celestial 
sphere with a number of evenly distributed spherical caps and 
then calculate ctj for each cap. The o-j obtained in this way 
for each cap can then be viewed as a measure of anisotropy 
in the direction of the center of that cap. Patching together the 
0" values for each cap, we obtain the indicator cr - cr{6, (p) 
defined over the whole celestial sphere, which measures the 
deviation from isotropy as a function of direction {9, <p). In this 
way, cr(6, (p) gives a scalar directional measure of anisotropy 
over the celestial sphere. 

Now, since a-{6, (p) is a discrete scalar function defined on 
S ^, we can further quantify the details of its anisotropy content, 
by expanding it in spherical harmonics in the form 

CO I 

(7(6, 4>)^YjY^ ^^'« 

f=0 m=-l 

and calculate the corresponding cr power spectrum 

in 

It then follows that, if a large-scale asymmetry is present in the 
original temperature distribution, it should significantly affect 
the cr-map on the corresponding angular scales (i.e. its lower 
multipoles). 

To recapitulate the multiple steps involved in defining our 
indicator cr, we itemize the procedure below: 

(i) Generate a set of evenly distributed overlapping spherical 
caps; 

(ii) For each cap, divide the set of pixels into a number of 
submaps A^jjjj, each containing an equal number of pixels 
within a temperature range; 

(iii) Calculate the PASH, •I'obs' ^'^^ ^^^ch submap; 

(iv) Average over all the ^Q\^g^ for each cap to obtain {<l>Q|^g); 

(v) Obtain the difference histogram T for each cap by subtract- 
ing Oexp - <'5obs>; 

(vi) Calculate the cr value, using Eq. (5), for each cap; 

(vii) Define the cr-map as the discrete function cr(0, 0) on 
where and (p are the coordinates of the center of each cap; 

viii) Calculate the multipoles De for the cr-map. 

From the above construction of cr{6, <p) and its power spec- 
trum, it is clear that there is no direct relation between the tem- 
perature multipoles and the o-(0, (p) multipoles. Although they 
are implicitly related, the relation is likely to be very compli- 
cated. In this work we take a more pragmatic approach by treat- 
ing our new anisotropy indicator as complementary to the usual 
temperature-based indicators, without attempting to establish a 
direct connection between them. 

In the next section we shall apply the indicator cr(0, 0) to 
both the first and three-year WMAP data. 
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Fig. 2. The full sky map of o-(6,(p) calculated from 12288 
evenly distributed 30° caps covering the ILC map. 

3. Application to WIVIAP data 

The WMAP measurements have produced high angular res- 
olution maps of the temperature field of the CMB. Here we 
considered both the first-year and the three-year WMAP data. 
For the first-year data we employed the full-sky foreground- 
cleaned Lagrange-ILC (LILC) map (Eriksen et al. 2004b), the 
TOH map (Tegmark et al. |2003] l. as well as the coadded WMAP 
map (Hinshaw et al. 200T). In all cases we chose HEALPix 
parameter A^gjjjg = 128 (Gorski et al. I2005I I. corresponding 
to a partition of the celestial sphere into 196608 pixels. For 
the three-year WMAP data we used the coadded map and the 
foreground-cleaned maps ILC (Hinshaw et al. l20061 l. OT (de 
OUveira-Costa & Tegmark'2006'), and PPG (Park et al. l2006] l. 

In all, we performed the cr-map analysis for seven CMB 
temperature maps: three using the first-year data and four with 
the full three-year data, all with different foreground cleaning 
algorithms. However, since the results are largely similar, we 
only present a detailed analysis of the three-year WMAP ILC 
map in the following. We also present the cr-maps for four of 
the remaining maps for comparison (the other two are similar 
to the depicted maps and were dropped to avoid repetition). 

In our calculations of the cr-map for these CMB maps, we 
subdivided the sphere into a set of A^caps - 12288 spherical 
caps of radius 30° co-centered with the same number of pix- 
els generated by HEALPix with A^gj^jg - 32. We subdivided 
each cap into N^[^ - 64 submaps, each consisting of the same 
number of pixels, and the number of bins chosen to produce the 
PASHs was taken to be A^^ins - ^^0. 

Figure |2] shows the Mollweide projection of the cr-map in 
galactic coordinates obtained from the ILC WMAP map. As 
can be seen, the cr-map shows a high-cr region suggestive of 
a dipole-like distribution. This result indicates that the distri- 
bution of hot and cold temperature pixels are not evenly dis- 
tributed in the celestial sphere. We note that this is different 
from a temperature dipole, where one half of the sky would, on 
average, be hotter than the other An anisotropy cr-dipole is not 
directly related to the temperature dipole, which was removed 
in all the CMB maps we have used, but is instead a result of 
large-angle anisotropics that are not readily apparent from a 
simple inspection of the temperature map alone. Particularly 
prominent is a very high-cr region concentrated on the south 



Fig. 3. The power spectrum of the cr-map for the ILC map, 
for { - 1, • ■ • , 10 (solid dots). Note the high values of the first 
three multipoles. For comparison, the mean expected power 
spectrum of the cr-map obtained by averaging over a set of 
1000 Monte Carlo-generated, statistically isotropic CMB maps 
(solid curve). Shown also are the corresponding cosmic vari- 
ance bounds for cr-maps (shaded area). 

eastern corner of the cr-map with a well-defined maximum at 
(b = 108°,/ = 222°). This is within a 16.6° separation from 
the directions recently indicated by Eriksen et al. (I2004al l and 
Land & Magueijo (I 2005a| l. 

To obtain more quantitative information about the observed 
anisotropy, we calculated the power spectrum of the cr-map for 
the ILC map, which is represented by solid dots in Fig. [3] To 
estimate the statistical significance of these cr-multipole val- 
ues, we compared the power spectrum of the cr-map, calcu- 
lated for the ILC map, with the averaged power spectrum of 
the cr-map obtained by averaging over a set of 1,000 Monte 
Carlo-generated statistically isotropic CMB maps. Considering 
the WMAP best-fitting angular power spectrum of the ACDM 
model (Spergel et al. I2006I I. each Monte Cai'lo CMB map is 
a stochastic realization of this spectrum obtained through ran- 
domized multipole components flf,„ generated within the cos- 
mic variance limits. 

The ( values used in this comparison lie in the range 1 < 
^ < 10, which are sufficient to allow the investigation of the 
large-scale angular correlations corresponding to angular sepa- 
rations > 18° in the cr-map. Moreover, we observe that Dc -0 
for ^ > 8. Averaging over the 1,000 spectra of the cr-maps 
obtained from the Monte Carlo CMB maps gives the expected 
power spectrum, which is depicted in Fig. [3] (continuous line) 
along with the corresponding cosmic variance! (shaded area). 

We note that the very high dipole value (corresponding to 
f = 1) found here is consistent with the cr-map (see Fig. |2]l 
showing a clear separation of the higher and lower values for 
cr into roughly two hemispherical regions. Using a test we 
have found this high value of the dipole to be statistically sig- 
nificant at 97% confidence level, compared to the expected 
isotropic power spectrum. Note also the high values of the 
quadrupole and octopole moments that were found to be sig- 

* We have checked numerically that the hi,„s of the cr-maps derived 
from the Monte Carlo CMB maps do follow Gaussian distributions, 
which justifies the use of cosmic variance bounds and the^- test. 
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nificant at 99% and over 99.9% (corresponding to a likelihood 
of one part in 10"^ in a statistically isotropic map), respectively. 
We also find that the likelihood of the first three multipoles of 
an isotropically-generated map to have such high values simul- 
taneously is very small being less than one part in 10^. Of the 
subsequent multipoles, only D4 and De fall outside the cosmic 
variance bounds. 

To check the robustness of our results we also calculated the 
power spectrum cr-map obtained from the WMAP PPG map. 
We found that the first three multipoles remain high relative to 
what would be expected from an isotropic case, at above 97% 
statistical significance. 

Given the anisotropic nature of the calculated cr-map, it 
is worth analyzing the shape of the anomalous multipoles in 
more detail. To this end, we depict in Fig. |4] the dipole, the 
quadrupole, and the octopole, as well as the full cr-map with 
these three components and the monopole removed. The anal- 
ysis of the dipole alone gives a direction towards (b = 130°, I = 
259°), which on its own does not agree very well with the axis 
of maximum asymmetry found by Eriksen et al. (I2004al) and 
also by Land & Magueijo ('2005a). However, the sum of the 
first three multipoles has two maxima, one close to the galactic 
center, at (b - 97°, I = 342°) and probably due to unremoved 
foregrounds, and the other in the direction {b - 1 12°, / - 217°), 
which is in much better agreement with the axis of maximum 
asymmetry indicated by Eriksen et al. ( I2004al) and Land & 
Magueijo (I2005al) . 

The quadrupole component has a very peculiar shape. It 
has two maxima and two minima lying symmetrically along 
a great circle. To determine this feature more precisely, we 
used the method proposed by de Oliveira-Costa et al. (I2003I I. 
where for each given multipole (in their case, of the tempera- 
ture map) the direction for which 'angular momentum disper- 
sion' is maximized is said to be a preferred direction. Applying 
this to the cr-map, we find that the angular momentum disper- 
sion for the quadrupole component is maximized in the direc- 
tion (b = 80°, / = 316°). 

The calculation of our anisotropy indicator requires not 
only the choice of a CMB map as input, but also the specifi- 
cation of a number of ingredients and parameters whose choice 
could in principle affect the outcome of our results. To test the 
robustness of our scheme, hence of our results, we studied the 
effects of changing various parameters employed in the calcu- 
lation of our indicator. We found that so long as the statistical 
noise (see Bernui and Villela i2006i) satisfies 

^^^^ < 0.05, (8) 

«#hir 

where A'jjjg is the number of submaps and n is the mean num- 
ber of pixels in each submap, then the cr-map and its power 
spectrum do not change appreciably as we change the number 
of submaps N^[^ from 2 to 64, the size of the spherical caps be- 
tween 15° and 30°, the number of spherical caps, A^caps (with 
values 768, 3072, and 12288, corresponding to HEALPix pa- 
rameter A^gj^jg - 8, 16, and 32, respectively) and the A^^ins 
ues in the range 180 to 400. 

Given that galactic foregrounds are likely to leave their 
most pronounced effects in the neighborhood of the galactic 



equator, one may suspect that the pronounced features, such as 
the high cr region in the south eastern corner of the cr-map (see 
Fig. (I2]i), may be an artifact of the foreground contamination 
processes. 

To show that this not the case, we first note that the lat- 
ter feature lies outside the exclusion region defined by the KpO 
mask, which is the severest cut-sky WMAP mask. To proceed, 
we also calculated the cr-maps for the other 6 CMB temper- 
ature maps. We studied the three-year (Hinshaw et al. I2006I I 
and one-year coadded (Hinshaw et al. 2003) maps, which are 
the combination of eight (from Ql trough W4) foreground- 
reduced sky-maps in the Q-, V-, and W-bands, and mini- 
mized but did not seek to eliminate the galactic contamina- 
tion. We also studied the foreground-cleaned TOH and OT 
maps given by Tegmark et al. (12003) and de Oliveira-Costa & 
Tegmark (2006), as well as the LILC (Eriksen et al. l2004bl l and 
PPG (Park et al. 2006) CMB maps. The corresponding cr-maps 
are shown in Fig.|5](the PPG and one-year coadded maps were 
dropped to avoid repetition). Together, these figures show that, 
despite differences in detail, the prominent high-cr region in 
the south eastern corner remains robust. 

These results indicate that the main anisotropic feature we 
have found is independent of the foreground cleaning processes 
employed. In addition we recently (Mota et al. 2007 ) performed 
a similar analysis of the WMAP one-year LILC data with both 
KpO and Kp2 masks applied and found very similar results. We 
also showed that this feature remains essentially the same in 
the Q, V and W bands (when analyzed after applying the Kp2 
mask) which would not be expected if it were an artifact of 
uncleaned galactic foreground. 

We should also note that apart from the high cr region in the 
south eastern corner, which lies outside the KpO mask, there is 
another weaker structure around the galactic center. Although 
this structure remains in some form in all the maps we have 
analyzed its intensity and other details vary considerably for 
different choices of the cleaning algorithms, as can be seen 
from Figs.|2]and|5] Considering that it lies in the most contam- 
inated region of the galactic plane, well inside the KpO mask, 
it is not clear whether this feature is not the result of the galac- 
tic foreground contamination. Furthermore its contributions to 
the excess dipole, quadrupole, and octopole is relatively minor 
compared to the high-cr region in the south-eastern corner. 

4. Conclusions 

We have proposed a new method of measuring directional de- 
viations from statistical isotropy in the CMB sky, in order to 
study the possible presence and nature of large-scale anisotropy 
in the WMAP data. 

Our anisotropy indicator has enabled us to construct a cr- 
map in order to search for evidence of large-angle anisotropics 
in the WMAP CMB temperature field. In particular we have 
found, with high statistical significance, a small region in the 
celestial sphere with very high values of cr, indicative of an axis 
of asymmetry that defines a direction close to the one reported 
recently (Eriksen et al. l2004al Land & Magueijo 2005al) . 

To obtain a more quantitative measure of this anisotropy, 
we studied the spherical harmonic expansion of the cr-map 
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Fig. 4. Depicted clockwise are the ILC cr-map's dipole, quadrupole, octopole, and the remaining cr-map from which all these 
components plus the monopole were removed. 
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Fig. 5. Depicted clockwise are the cr-maps: LILC, three-year coadded map (showing the effects of the unfiltered galactic con- 
tamination), OT and TOH WMAP maps. The one-year coadded and the PPG maps are not shown as they are very similar to the 
three-year coadded and the ILC maps, respectively. Note the high-cr region in the south eastern corner of the cr-maps, which 
remains robust for all maps. 



generated from the ILC map and found that the correspond- 
ing dipole component is significantly higher than would be ex- 
pected, as compared with an average of a set of cr-maps ob- 
tained from 1000 Monte Carlo CMB maps generated under the 
statistical isotropy hypothesis. Moreover, we have found that 
the combination of the dipole, quadrupole, and octopole com- 



ponents of the cr-map has a strongly preferred direction close 
to the small region of the maximum value of cr. 

Using a test we found that this high value of the dipole 
is statistically significant at a 97% confidence level (CL), com- 
pared to the expected isotropic case. We also found high values 
of the quadrupole and octopole moments, which were found 
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to be significant at 99% CL and over 99.9% CL (correspond- 
ing to a likelihood of one part in lO'* for statistically isotropic 
generated maps) respectively. The likelihood of the first three 
multipoles of an isotropically generated map to have such high 
values simultaneously is very small at less than one part in 10^. 

We have shown that the results reported here are robust by 
showing that the cr-map does not significantly change with dif- 
ferent choices of parameters, so long as the statistical noise is 
kept under control. We have also studied the effects of different 
foreground-cleaning algorithms, or absence thereof, by consid- 
ering the LILC (first-year), TOH (first-year), OT (three-years), 
PPG (three-year), and the coadded (first and three-year) CMB 
maps in addition to WMAP three-year ILC map. 

We have found again that the corresponding cr-maps re- 
main qualitatively unchanged, in that the high-cr region in the 
south eastern corner of the cr-map remains essentially invariant 
for all the maps considered here, while similar statistical sig- 
nificances for the most important (low-^) multipoles D( hold 
for all CMB maps considered. This provide evidence that our 
results regarding the large-angle anisotropy in the CMB maps 
remain unchanged no matter whether the first or three-year data 
are used. 

This robustness demonstrates that our indicator is well- 
suited to the study of anisotropics in the CMB data. 
Furthermore, although derived from the temperature maps, it 
reveals features that are not apparent in the analysis of the 
temperature maps alone. Our scheme for studying directional 
anisotropics is complementary to the existing approaches in 
the literature, in the sense that it reveals different aspects of 
the same data set. However, the relationship between them is 
clearly non-trivial and requires further investigation, which is 
beyond the scope of the present work. 

Regarding the origin of such large-angle anisotropy, a num- 
ber of suggestions have been put forward. Briefly, they can 
arise either from a subtle form of unremoved foreground con- 
tamination (in which case the cr-map might indicate where in 
the sky this contamination is most intense) or from the universe 
being genuinely anisotropic on large angular scales. This pos- 
sibility is particularly interesting, as it would have potentially 
important consequences for the standard inflationary picture, 
which predicts statistically isotropic CMB temperature fluctu- 
ation patterns. 

Among this type of suggestions, it has been proposed that 
the preferred direction could be due to the universe possessing 
a non-trivial topology (see, e.g., Copi et al. 120051 Spergel et 
al. l2003t for more details on cosmic topology see, e.g., the re- 
view articles Lachieze-Rey & Luminet |l995l Starkman |1998l 
Levin l2002t Rebougas & Gomero l2004t Lehoucq et al. [19961 
Roukema & EdgeHlHl Corni sh et a l.[l998l Bo nd eta l.|2000il 
l2000b; Fagundes & Gausm ann|T999l Uzan et al. [T999l Gomero 
et al. l2000ll200Tall2001bll200Tcl 120021 Mota et al. l2003ll2004l 
Luminet et al. l2003l Aurich et al. '2005a 2005b; Hipolito- 
Ricaldi & Gomero 12005. Weeks et al. ,2003a. Weeks ,2003b ). 
If topology is indeed the origin, our anisotropic indicator is 
promising for distinguishing between different topologies. On 
the other hand, the extent to which this is true for the cr- 
indicator is unclear and deserves further investigation. 



Whatever the origin of these large-angle anomalies may be, 
the robustness of our anisotropy indicator seems to be suffi- 
ciently sensitive to reliably map them, which in turn could fa- 
cilitate the task of explaining their origin. 
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